Can scalar leptoquarks explain the f Ds puzzle? 
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Motivated by the disagreement between experimental and lattice QCD results on the D s decay 
constant we systematically reinvestigate role of leptoquarks in charm meson decays. We consider 
scalar leptoquarks that transform as a weak interaction triplet, doublet, or singlet in a model 
independent approach, and also argue that in a particular SU(5) GUT model these leptoquark 
states, contained in the 45-dimensional Higgs representation, could be safe against proton decay 
bounds. Using the current experimental measurements in r, kaon and charm sectors, we find that 
scalar leptoquarks cannot naturally explain the D s — > [iv and D s — > tv decay widths simultaneously. 
While any contributions of the triplet leptoquarks are already excluded, the singlets could only 
contribute significantly to the D 3 — > tv width. Finally, a moderate improvement of the experimental 
upper bound on the D° — > decay width could exclude the doublet contribution to the D s — > 

fiv, while present experimental data limits its mass to be below 1.4 TeV. Possible new signatures 
at present and near future experiments are also briefly discussed. 

PACS numbers: 13.20.Fc,12.10.Dm,12.15.Ff 



I. INTRODUCTION 

Leptoquark states are expected to exist in various ex- 
tensions of the Standard Model (SM). They were first in- 
troduced in the early grand unification theories (GUTs) 
in the seventies [1, 2]. Scalar leptoquarks are expected to 
exist at TeV scale in extended technicolor models as well 
as in models of quark and lepton compositcness. Scalar 
quarks in supcrsymmetric models with R-parity viola- 
tion (RPV) may also have leptoquark-type Yukawa cou- 
plings [32]. 

Recently, discrepancies between the experimental mea- 
surements of leptonic decay modes of D s mesons [3-6] 
and the lattice results for the relevant fu 3 decay con- 
stant [7-9] have stimulated many analyses. One intrigu- 
ing indication is that the central measured and predicted 
values for the fo s differ by more than 10 % with a com- 
bined significance of D s — > tv and D s [w channels 
of roughly 2.3 c [10], while the corresponding values for 
f D are in perfect agreement. In Rcf. [11] the idea of 
scalar leptoquarks has been revived to explain the miss- 
ing decay widths. Some implications of this suggestion 
have been further explored using semileptonic [12, 13] 
and rare charm decays [14]. 

Generally, leptoquarks which also couple to diquarks 
mediate fast proton decay and are therefore required to 
be much above the electroweak scale [15], making them 
uninteresting for other low energy phenomena. "Gen- 
uine" leptoquarks on the other hand, couple only to pairs 
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of quarks and leptons, and may thus be inert with respect 
to proton decay. In such cases, proton decay bounds 
would not apply and leptoquarks may produce signatures 
in other low-energy phenomena. In this article we set 
out to study whether scalar leptoquarks can naturally 
account for the fn s puzzle and at the same time comply 
with all other measured flavor observables. 

We consider all possible renormalizable leptoquark in- 
teractions with SM matter fields consistent with the SM 
gauge symmetry. One can construct such dimension- 
four operators using leptoquarks which are either sin- 
glets, doublets or triplets under the SU(2) L . If we fur- 
thermore require that such leptoquarks contribute to lep- 
tonic decays of charged mesons at tree level, we are left 
with three possible representation assignments for the 
SU(3) C x SU(2) L x U(1) Y gauge groups: (3,3,-1/3), 
(3,2,-7/6) and (3,1,-1/3). Only the weak doublet 
leptoquark is "genuine" in the above sense. However, 
using a concrete SU(5) GUT model where the relevant 
leptoquarks are embedded into the 45-dimensional Higgs 
representation (45 h), we demonstrate how leptoquark 
couplings to matter can arise and in particular, how the 
dangerous couplings to diquarks - both direct and indi- 
rect [15] - can be avoided. 



II. GENERAL CONSIDERATIONS 

In our analysis we will assume the mass cigenstates 
within a leptoquark weak multiplct to be nearly degen- 
erate. While large mass splittings within a weak mul- 
tiplet may be considered unnatural, more importantly, 
they are also tightly constrained by electroweak preci- 
sion observable T [16]. Consequently, one gcncrically 
gets correlations between semileptonic charged currents 
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and (lepton) flavor violating neutral currents, which rep- 
resent important constraints on any leptoquark scenario 
trying to resolve the D s leptonic widths puzzle. Also, 
we focus on observables mediated by the relevant lepto- 
quark couplings at tree level since these already involve 
processes forbidden in the SM at tree level, i.e., flavor 
changing neutral currents (FCNCs) and lepton flavor vi- 
olation (LFV) processes. Finally, since the present fo s 
deviation is of mild significance, we require all the mea- 
sured constraints to be satisfied within one standard de- 
viation (at 68% C.L.) except upper bounds, for which 
we use published 90% C.L. limits. We consider a lepto- 
quark explanation of the fo a discrepancy as natural, if 
both D s and D leptonic decay widths can be obtained 
close to their measured central values. 

After the clcctrowcak (EW) symmetry breaking, 
quarks and leptons acquire their masses from their re- 
spective Yukawa interactions. These induce in the phys- 
ical (mass eigen-)basis a CKM and PMNS rotations be- 
tween the upper and the lower components of the fcrmion 
doublets. Consequently, it is impossible to completely 
isolate leptoquark mediated charged current interactions 
to a particular quark or lepton generation in the left- 
handed sector irrespective of the initial form of the lepto- 
quark couplings to SM matter fields, unless there is some 
special alignment with the right-handed quark sector. To 
see this, we denote as J''' a 3 x 3 arbitrary Yukawa ma- 
trix in the weak basis, and write down flavor structure of 
interaction of the quark and lepton doublet parts 

QfXi e =(^ d^)X^ = {w q dT q )(U^Xy e , (la) 
X' qe Lf =X'«V7 e^) T = (X'E) qt (i/ t e t ) T , (lb) 

where fields with w superscript are in the weak basis, 
whereas d 1 = VcKud and v' = Vpmns v - The uni- 
tary matrices U, D, E, and N rotate the fields from mass 
to weak basis and are unphysical per se, so we absorb 
them in redefinition of the couplings (e.g. Ylq = W X 
on the quark and Y' L q = X'E on the lepton side) and 
consider them as free parameters. In what follows, we 
will use the convention where all the remaining rotations 
are assigned to down-type quark (Vckm = U^D) and 
neutrino (Vpmns = E^N) sectors, and the quark mass- 
eigenstates are defined as 

(Qi,Q 2 ,Q 3 ) - s c , 6 *) , K s' V) = (d s b) v£ KM . 

It is obvious in this notation that, even if the Ylq matrix 
had all rows, except for the g-th one, set to zero, which 
would correspond to leptoquark coupling only to u q , one 
would still get non-zero couplings to all three left-handed 
down-quarks. Same rationale holds true for the lepton 
sector due to Vpmns, but since the neutrino flavors are 
not tagged in present experiments, the respective decay 
widths are summed over all neutrino flavors. Whenever a 
mass-eigenstate antineutrino vi is produced in a reaction, 
its amplitude includes, according to Eq. (1), a factor of 
^j^LQ^PMNS f° r leptoquark interaction, or Vp MNS if 



the neutrino was produced in W tv vertex. In any case, 
when one sums the rates for all neutrino species 

l-4i| 2 ~ VpMNsVpMNS — & l ! 

j=l,2,3 i=l,2,3 

it becomes evident that in the summed rate all the neu- 
trino indices are replaced by the lepton flavors. This is 
equivalent to the absence of mixing in the lepton dou- 
blets. 

The above considerations are more general and similar 
in spirit to the ones recently discussed in Ref. [17] for the 
case of K—K and D — D mixing. In fact, any new physics 
coupling to SM fcrmionic weak doublets exhibits similar 
kind of correlations, and contributions to charged cur- 
rent transitions cannot be isolated to a particular quark 
generation. 

Another important particularity of the fu 3 puzzle is 
that it affects a Cabibbo favored c — > s transition. Con- 
sequently, the hierarchy of correlations with other pro- 
cesses is largely determined by the CKM mixing hierar- 
chy. In particular, the mixing of the third generation 
with the first two is much smaller than the mixing of 
the first two generations among each other. Therefore, 
for our purposes, it is often a good approximation, to 
completely neglect effects of the third generation in the 
quark sector. Then we can parameterize a generic lepto- 
quark coupling in the weak basis using a common (real) 
pref actor and a rotation angle 

Y lq = 2/lq (sin 0, cos <j)). 

In addition, the only CKM rotation is due to the Cabibbo 
angle and there is no SM CP violating phase (d! — 
cos 9 c d+siii 9 c s, s' = — sin 9 c d+cos 9 c s and V us = —V c d = 
sin# c = 0.225). The absence of SM phases is not critical 
for our purposes, since we only consider CP conserving 
quantities, and since the relevant SM amplitudes in our 
considered processes have approximately the same weak 
phase even in the full three generation case. The lepto- 
quark couplings themselves, however, could in principle 
have arbitrary new phases. These could be important in 
processes with two or more interfering amplitudes con- 
tributing, at least one of those being due to the lepto- 
quarks. We deal with this possibility on a case by case 
basis. Finally, in all the scenarios considered we have 
checked explicitly that the two generation approximation 
is valid by performing numerical leptoquark parameter 
scans including the full CKM structure and a full set 
of possible leptoquark couplings with arbitrary phases. 
In this case, (semi) leptonic B decays B — > tv and espe- 
cially B — > Dtv can be used to put additional constraints 
on the leptoquark parameters relevant to the D s — > tv 
width. Numerically however, these constraints turn out 
not be competitive with the others due to presently lim- 
ited experimental precision. 
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III. TRIPLET LEPTOQUARK (3,3,-1/3) 

The triplet leptoquark can in principle couple to di- 
quarks and thus destabilize the proton, so one has to 
check in an underlying model if that is indeed the case. 
The allowed leptoquark interaction Lagrangian consists 
of a single term 



A3 Lj + h.c. . 



(2) 



where Q° = -Q T C^ 1 , C = i-/ 2 ^ and r are the Pauli 
matrices. The 3x3 coupling matrix Y 3 is arbitrary in the 
bottom-up approach. On the other hand, its entries may 
be related to other parameters in an UV embedding of the 
effective theory. In the concrete SU(5) model analyzed 
in the Appendix A, the above couplings are due to the 
contraction of 10 and 5 with 45* H - also responsible for 
giving masses to the down quarks and charged leptons. A 
different contraction of 10 and 10 with 45# couples the 
triplet to diquarks. The latter term can be consistently 
set to zero in the supcrsymmetric version of the model, 
thus sufficiently suppressing proton decay. 

As already mentioned in Ref. [11], the triplet lepto- 
quarks cannot by themselves account for deviations in 
D s — * tv for both t and (i in the final state due to con- 
straints coming from LFV tau decays, such as r — > rf' V 
and r — > (f)fi. Numerically, the r — > ijfi decay turns out 
to be most constraining. The triplet leptoquark contri- 
bution can be written as 



r(3) = 

x t— >rju 



E,. T><3l"T>qM* fq 
q=u,d,s Xg I 3 1 3 Jti 



5127rmi 3 



-m" 



m 7 



-1 3/2 



where we have neglected the muon mass. The weight 
X q = 1 is for q = u and 2 for q = d,s comes from 
an additional \[2 factor in the interaction terms with 
^3(^3 = il) states. Decay constants of 77 meson /' are 
defined as in [18]. As explained in the previous section, 
the couplings Y 3 contain an additional Vckm rotation 
for the down-type quarks 



ql 



{VckmYzY 



q = u,c, t, 
q = d,s, b. 



Thus, the upper bound on r 



rjfi decay width 



directly constrains the products Y 3 ST Y 3 ^* , Y 3 T Y 3 ^* 
and Y 3 UT Y 3 utl *(=(cos9 c Y 3 dT + sin 6> c F 3 sr )(cos 6 C Y^* + 



smO c Y 3 st1 *) in the two-generations approximation). 



On 



the other hand, the relative contributions to the D s lep- 
tonic widths (see Eq. (3)) in the two generations approxi- 
mation are Y 3 Sfl (Y 3 ^* —ta,n9 c Y 3 dt1 *) for the muon channel 
and Y 3 ST (Y 3 ST * - tan 6 c Y 3 dT *) for the tau channel. Both 
cannot be sizable and at the same time agree with the 
bounds that come from the r — > r\\x decay width [33]. In 
scenario of triplet leptoquarks therefore, one of the mea- 
sured leptonic channels D s — > Iv would necessarily have 
to be a measurement artifact. We will consider both pos- 
sibilities separately. 
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If the leptoquarks have sizable coupling Y 3 ST (implying 
by the r — > 77/4 decay width) we can obtain a 
non-zero contribution to the D s — ► tv decay width due to 
the interference term between the SM and the leptoquark 
amplitude in 



p( 3 ) _ r SM 

1 D S ^TV ~ 1 D S ^TV 



where the SM width 



sm _G 2 F m 2 T \V cs \ 2 f 2 D m D 



8tt 



m D s 



and the relative contribution of the triplet leptoquarks 
reads 



Yi 



r 3 



V cs m 2 A 



(3) 



An important observation is that the leptoquarks in this 
scenario contribute to the same effective operator as the 
SM and thus exhibit the same helicity suppression. In 
the two generations approximation, the relative triplet 
contribution simplifies to ^3 = (y 3 ) 2 cos <j>sm </>(tan 9 C + 
cot 4>) I to^ 3 . Reproducing the measured branching ratio 
Br(D s -» tv) = 0.0561(44) [10] while using the most 
precise lattice input fo s — 241(3) MeV [7] would re- 



quire 



0.002 GeV" 



On the other hand the Y 3 T 



coupling of leptoquarks is constrained by precise mea- 
surement of the lepton flavor universality ratio Tr^/ T = 
Br(r -> TTv)/Br{TT -» fiv) = 0.1092(7) [19, 20]. Lepto- 
quarks contribute to semileptonic tau decays in the form 



r (3) 



where 



-.SM 



1 + 



1 

4x/2Gi 



V ud m\ 



(4) 



-iSM 



G 2 F ml\V nd \ 2 f 2 
16tt 



m 7 



In the two generations approximation the term in square 
brackets can be written as (yj) 2 sin </>cos </>(— tan C + 
tan (f>) /m\ a . To exactly satisfy both the tt^/t value and 
explain leptonic D s — > tv excess we need either 

(a) tan</> w tan6* c , i.e., leptoquarks couple only to s 



but not to d quark (Yj-. 



dr 



0), 



(b) sin</> i=a 0, i.e., leptoquarks couple only to c but not 
to u quark (Y 3 UT = 0). 

However, in the limit (a) one must have a sizable cou- 
pling Y 3 UT due to CKM rotation which results in rela- 
tive contribution of size S T to the Cabibbo suppressed 
semileptonic tau decays r — > Kv (of the form (4) with 
appropriate flavor replacement d — > s and tt — ► K). 
These are measured in agreement with the SM at the 



4 



3% level [34] (in particular the ratio -K^/t = Br(r — > 
Kv)/Br(K -> j*z/) = 0.0109(4) [19, 20]) [35]. 

In the other limit, (b), one must have sizable coupling 
and thus gets a relative contribution scaling as S T 
to the Z? — > tv decay width. Currently only an upper 
bound exists for this channel Br(D — > tv) < 1.2 x 10~ 3 
at 90 % C.L. [4]. Even more importantly, one gets a non- 
vanishing contribution to the rare K + — ► 7r + z/P decay. 
Since the triplet leptoquark contributes with the same 
effective operator as the SM, its contribution can be ob- 
tained by simply replacing the A t A t product in the mas- 
ter formula of [21] with 



XtX t 



XtX t 



GFOt em Sin 0\ylTi^ 



This process is measured to have Br(K + 



~vv) 



(l7.3t{o%) x 10_11 I 22 > 23 ]- :t constraints the 
sum of leptoquark coupling combinations Y^Y^* = 
(y|) 2 cos C cos 0(tan — tan# c )(l + tan C tan (f>) and fixes 
very accurately tan</> = tan# c or tan0 = — cot9 c . The 
constraint applies to all lepton flavors since it is in- 
clusive with respect to neutrino flavor [36]. The com- 
bined impact of all these constraints on the triplet lep- 
toquark contribution to D s — > tv is shown on the first 
plot in Fig. 1. One observes that the combination of 
the strong bounds coming from K + — > tt + vv combined 
with Br(r — > Kv) / Br(K — > \xv) completely excludes the 
triplet lcptoquarks from explaining the D s — > tv excess. 

In the opposite scenario where the leptoquarks couple 
to muons the situation is similar to the tau case with 
two differences: (1) the D — > /iv decay width has already 
been measured and the Br(D — > \iv) = 3.8(4) x 10 -4 [4] 
agrees perfectly with the SM prediction using the most 
precise lattice QCD value of f D = 208(4) MeV [7]; (2) an 
additional constraint comes from the FCNC decay Kl — > 
as it receives contributions from leptoquarks of the 

form 



r 



(3) 
K L 



fKm K 7 



1 



64tt 



l A 3 



The requirement that such leptoquark contributions do 
not exceed presently measured Br{K^ — > /i + // _ ) = 
6.84(11) x 10 -9 [19] produces a bound equivalent to the 
existing one coming from K + — > ir + vv. Combining these 
two additional constraints with the rest also clearly dis- 
favors a triplet leptoquark explanation of the D s — * [iv 
excess, as shown on the bottom plot in Fig. 1. 



IV. DOUBLET LEPTOQUARK (3,2, 7/6) 
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Figure 1: Combined bounds on the triplet leptoquark param- 
eters in the two-generation limit in the tau (upper plot) and 
muon (lower plot) sectors. All bands represent 68% C.L. ex- 
clusion intervals, except the upper bound on D — > tv which 
is taken at 90% C.L.. The K + -> ty + vv and Kl -> mV" 
constraints can only be satisfied on the two horizontal dashed 
lines. Within the green bands, the D s — + tv excess can be 
accounted for. 



In the particular SU(5) model, the term proportional to 
Y 2 r stems from the contraction of 10 and 5 with 45^ 
while the Y 2 l term is due to 10 and 10 being contracted 
with 45 h- 



The doublet leptoquarks are innocuous as far as proton 
decay is concerned. The allowed dimension four interac- 
tions in this case are 

£2 = Y 2 l i Q, it 2 A* ej + Y 2 % UiAlLj + h.c. . (5) 



In this scenario the same states couple left-handed 
quarks to right-handed leptons and vice versa. Con- 
sequently, only the product of both couplings can con- 
tribute to the D s leptonic widths through the interfer- 
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ence with the SM 



-pSM 

1 D s ^eu 



4x/2Gj 



(6) 



I 2R X 2L 



fit = _ u, 

2 m e (m c + m s ) V* s m\ 



Again, the couplings of left-handed down- and up- type 
quarks are misaligned 

Y 2 i = (v£ KM Y 2L y e forq = d,s,b. 

Note that the doublet leptoquark contribution exhibits 
no helicity suppression. Thus, explaining both muon and 
tau leptonic partial widths of D s requires vastly different 
leptoquark couplings. On the other hand, now one also 
has to take into account the strict bound coming from the 
decay D° — » In the doublet leptoquark model, 

this mode receives potential contributions from several 
coupling combinations 



fp™p 

' 5127T 



Am 2 , 



rn 



D 



4m 2 \ 

— # \A + B 



(7) 



where A and B contain the couplings of doublet lepto- 
quarks 



m% (m c + m u ) 2 1 2R 2L 



B 



1 



\ I 2R 1 2L 



m c + m„ 



2m M 
m D 



V2R J 2 



Ufl 



1 2R 1 2L\ 



v cfj,* v ufj,\ 
1 2R 1 2L I 



Two combinations involve Y 2 £ which is in conjunc- 
tion with Y 2 £ and Y 2 [ constrained through preci- 
sion kaon and tau lepton flavor universality tests 
similarly as in the triplet scenario. In addition, 
this coupling does not contribute to the D s — > [w 
width (6). The remaining two combinations can be 
rewritten by using the Cabibbo rotation in terms of 



1 2R 1 2L 



Y^(cos9 c Y 2 d r+sm8 c Y 2 T) and Y%Y 2 T = 



(cos 6 C Y 2 £ - sin 6 C Y$) (cos C Y 2 £ 



rdfJ.* 



sin 6 C Y 2 £* 



The 



D — > fiv width measurement constrains directly the 
size of Y 2 j* . In absence of this coupling, D° — > 
would receive dominant contributions from just two non- 
interfering leptoquark amplitudes 



/?? ■ 



B 



sin 2 6 r 



- m u ) 
m D 



1 2R 1 2L bmt7 c 



~ I 2R I 2L 



2m 



; V 2L ) 



cos 9, 



where A can be related to D s — > /iv decay width con- 
tribution (it is proportional to sin 2 6> c cos 2 C \6 2 | 2 ). The 
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Figure 2: Combined D° — -> (i + \i~ and D — > fj,u bounds on the 
doublet leptoquark parameters in the two-generation limit in 
muon sector as explained in the text. D — ^ [iv band represent 
68% C.L. exclusion interval, while the upper bound on D° — > 
H + fi~ is taken at 90% C.L.. The Kl — * constraint can 

only be satisfied on the two horizontal dashed lines. Within 
the green band, the D s — + tv excess can be accounted for. 



term proportional to A alone yields Br(D° — > w 
8.3 x 1CP 7 for the central value of D s — > fiv decay 
width. Recently, an improved experimental limit of 
Br(D° -> < 4.3 x 1CT 7 at 90% C.L. was put 

forward by CDF [24]. It is evident that this introduces 
some tension between explaining the D s excess and not 
spoiling the agreement in the D case. Due to the mod- 
erate significance of the D s discrepancy, this tension is 
not yet conclusive as can be seen on Fig. 2, where we 
plot the combined constraints in the <f>-y 2 plane. Fig. 2 
is generated in the following way. We first parameter- 



ize y 2 "£ 



V2L C0S 1 



and Y- 



d/i 



YZ = while Y% 



2L 



V2L Sm < 



We then set 



y 2 R. Finally, we vary y^ L and y^ R 
while keeping the product y 2 — \J y 2 ^y 2 ji fixed and use 
the best fit value to determine the allowed region. We 
include the Bt(Kl — > fi + fi~) constraint which is also 
relevant in this case 



r (2) 



1 2L 1 2L 



fx m K 



1 



4m?, 



0-3) 



2567rrn^ 2 (m d + m s y 



While both central values for D and D s leptonic widths 
clearly cannot be reproduced by the doublet leptoquark 
contribution, a future improvement of the bound on 
D° — ► is clearly sought after to reach a definite 

conclusion on this scenario. 

Opposed to the triplet leptoquark case, the verdict on 
the D s — > tv contribution of the doublet leptoquark 
is still far from conclusive. Firstly, because the corre- 
sponding D leptonic mode has not been measured. Sec- 
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ondly, because there are presently no strong experimental 
bounds on FCNCs in the up quark sector involving only 
tau leptons or only neutrinos (doublet leptoquark does 
not contribute to s — > dvv transitions) . We note in pass- 
ing that, provided the doublet leptoquarks are to explain 
both tau and muon final state excesses, there is an impor- 
tant bound coming from the aforementioned r — > f/'V 
decays. In this scenario they constrain the following com- 
bination of parameters b^ FV = |^2L ^2L M l/ m A 2 appear- 
ing in 



p(2a) 



I 2L I 2L 



m 



I28irmh 



where we have only considered the s-quark contribution 
and have neglected the muon mass. On the other hand, 
explanation of D s leptonic excesses requires nonzero val- 
ues for $2 ' T - Finally, in order not to spoil perturbativc 
treatment of the couplings, none of the couplings should 
exceed a value of roughly \Y 2 t3 L R \ < \fH. Then, one can 
combine the above inequalities to yield a robust upper 
bound on the doublet leptoquark mass: 



m A2 < 



Taking the present Br(r -> < 6.5 xlO -8 at 90 % C.L. 
bound [19] and the central values for the two leptonic 
decay widths, one obtains a value of roughly 1.4 TeV, 
which is certainly within the LHC reach [25]. 



V. SINGLET LEPTOQUARK (3, 1, -1/3) 



where Y qt arc defined as in the triplet leptoquark sce- 
nario. The second term in Eq. (8) can come from the 
SU(5) embedding without causing any conflict with the 
bounds on proton decay lifetime even if the leptoquark is 
very light, whereas the presence of the first term would 
require some fine tuning in order for the leptoquark not 
to couple to diquarks (see Appendix A) . Note that if the 
first term is absent, then the singlet leptoquark cannot 
contribute to the D s leptonic decay width. If the second 
term is absent the analysis is analogous to the triplet 
leptoquark scenario, with the exception that the singlet 
does not contribute to Kl — * _ - Such is for example 
the case of the RPV minimal supcrsymmctric SM, where 
the interaction term of a right-handed down squark to 
quark and lepton doublets is present and corresponds to 
the first term in (8), while the second term is absent. 

From the triplet scenario we know that K + — > -k + vv 
forces the Y( L couplings to be diagonal in the down-type 
quark basis and in particular ss 0. Also relevant is 
the constraint from the lepton flavor universality ratio 
Br(r — > Kv)/Br(K — > fiv) which receives relative lepto- 
quark contributions of the form (9) with suitable flavor 
replacement (c — > u), while the tau semileptonic width is 
given as 



F 



(i) 



*Kv 



1 

4V2G f 



V... 



m 7 



, I 1R I IL 



V* s (m u + m s ) 



This state was originally proposed to explain the D s 
leptonic width puzzle in Rcf. [11]. On the other hand, 
singlet leptoquarks are notorious for their mediation of 
proton decay. However, as in the case of the leptoquark 
triplet, one can demonstrate in concrete SU(5) embed- 
ding (see Appendix A) that the dangerous couplings do 
not necessarily appear. Similarly to triplets, weak sin- 
glets can also couple to pairs of SM matter weak dou- 
blets. However, now also couplings to pairs of singlets 
are possible resulting in the dimension four interaction 
Lagrangian with two terms 

A = y/£Qfir 2 A^ + Y^Ale, + h.c. . (8) 

In this generic effective theory description clear correla- 
tions among different charged and neutral current flavor 
observables, present in the triplet case, are somewhat 
diluted by the presence of the second interaction term 
which modifies the singlet leptoquark contribution to the 
D s leptonic width 



,(i) 



1 + 



ycT*ysT 



AV2G F m 2 Ai 



V* s m T (m c + m s ) 



(9) 



Remaining constraint is the rare decay D° — > 
which in this case is of the form (7) with A and B 



''D 



B = 



m 



Ai 



m D 



m c + m. 
2m 



f v ufj, v cfi* 
\ 2 1R I 1L 



v cfj,* v ufj,\ 
I 1R I 1L I 



Y c t**Y uf *) 

\R 1L ) 



m D 



The remaining relevant free parameters can correspond- 
ingly be chosen as an overall coupling magnitude 5 and 
two angles (4>, u>), defined through Y*£ — y± sin cj, = 
Ui cos u cos <j> and Y^ = yf cos u> sin </>. The value of is 
bounded from above by the condition of perturbativity 
(■tji < \/47r). Together with existing direct experimen- 
tal searches for second generation leptoquarks [26, 27] 
this gives an additional constraint on the possible size of 
the leptoquark contributions to the D s leptonic width. 
By performing a numerical fit of (y^,<jj,<p) to these con- 
straints we obtain the result, that the experimental value 
for Br(D s — > fiv) cannot be reproduced within one stan- 
dard deviation without violating any of the other con- 
straints, thus excluding the singlet leptoquark as a nat- 
ural explanation of the D s — > \w puzzle. 
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As in the doublet case, the lack of experimental infor- 
mation on up-quark FCNCs involving only tau leptons 
leaves the verdict on the singlet leptoquark contribution 
to the D s — > tv decay width open. What is certain is that 
due to the K + — > n + vv constraint any such contribution 
has to be aligned with the down-type quark Yukawas such 
that Yf r i=a can be ensured. In addition, the second 
term in Eq. (8) needs to be present and sizable to avoid 
the bounds coming from Br(K — > /j,v)/Br(T — > if;/). 

VI. CONCLUSIONS 

Scalar leptoquarks cannot naturally explain both en- 
hanced D s — ► £f decay widths due to existing constraints 
coming from precision kaon, tau, and D meson observ- 
ables. The triplet leptoquark is excluded from contribut- 
ing to any of the widths. Sizable contributions due to 
single right-handed down squark exchange in RPV su- 
persymmetric models arc also excluded, while a generic 
leptoquark singlet is definitely excluded only from ex- 
plaining the D s — > fiv width. The doublet contribu- 
tion to this process is still technically allowed, while an 



improvement in the search for D° — > could very 

soon also completely exclude it. For the D s — > tv only 
the triplet (and RPV) explanations are already excluded, 
while the possible doublet explanation of both widths re- 
quires its mass to lie below 1.4 TeV and will certainly 
also be probed with direct leptoquark production at the 
LHC. Possible future signatures of a scenario where lep- 
toquarks are responsible for the D s — > tv width could 
also be Br(J/ip — > t + t~) at the level of 10 -11 , probably 
beyond the reach of BESIII [28], and also Br(t -» cr+r") 
at the level of 10~ 5 , close to the limiting sensitivity of the 
LHC [29]. 
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Appendix A: SU(5) EMBEDDING 

We now demonstrate i) how natural it is for the weak 
triplet, doublet and singlet leptoquark interaction terms 
to arise in renormalizable SU(5) model, and ii) how plau- 
sible it is for them to be light enough to play role in flavor 
physics phenomena. 

In SU(5), an zth (i = 1,2,3) generation of the SM 
matter fields comprises Wi(— (1, 1, l)i © (3, 1, — 2/3)j © 
(3,2,1/6), - (ep.uf.Qi)) and 5;(= (1,2,-1/2), 
(3, 1, 1/3), = (Li,df)), where Q l = (u; di) T and L; = 
{yi ei) T . The up quark (down quark and charged lep- 
ton) masses originate from the contraction of 10^ and 10^ 
(5j) with 5- and/or 45-dimcnsional Higgs representation. 
(Observe that 10x10 = 5©45©50 and 10x5 = 5©45.) 
Only these two representations contain component that 
is both electrically neutral and an SU(3) C singlet that 
can thus obtain phcnomcnologically allowed vacuum ex- 
pectation value (VEV). Actually, both are needed in a 
realistic renormalizable setting on purely phenomenolog- 
ical grounds. 

The most general renormalizable set of Yukawa cou- 
pling contractions with 5h and 45# is 

V = Yihof5 aj 5* H0 + Yi j e a0J seloflO] s 5h 

+ Y 4 %wf5 5j 45* H 5 a/3 + YHe a01 selof lof 45%, 

where Greek indices are contracted in the SU(5) space. 
Relevant fermion mass matrices are 

M D = (y£v* b + Tfl, v* 45 ) /V2, (Ala) 
M E = (Y 5 -*v* 5 - 6Y 45 .v* 45 ) /yfi, (Alb) 
M v = [4(F 5 T + Y 5 )v 5 - 8(Yg - Y 45 )v 45 ] /yfi, (Ale) 

where (5 5 H ) = v 5 /V2, (46&> = (45&) = (45f 3 > = 
v 45 /V2 and \v 5 \ 2 + \v 45 \ 2 = v 2 (v = 247 GcV). Y 5 , , Y 45 , , 
I5 and I45 are arbitrary 3x3 Yukawa matrices. 

If only 5-dimensional (45-dimcnsional) Higgs represen- 
tation were present one would have = (— 3)M D . A 
scenario with only one Higgs representation would hence 
yield m T /m(, = m^/m s = m e /md at the GUT scale, 
which is in conflict with what is inferred from experi- 
mental observations. This is why both 5 H and 45 H arc 
needed at renormalizable level. (Note, since m M /m s > 1 
whereas m T /m,b « 1 at the GUT scale this would suggest 
that the the Y 22 » entry is enhanced compared to other 
entries of Y 45 , [30].) 

Conveniently enough, the 45-dimcnsional Higgs rep- 
resentation 45h(— (Ai, A 2 , A 3 , A 4 , A 5 , A 6 , A7) = 



(8, 2, 1/2) ® (6, 1, -1/3) ® (3, 3,-1/3) © (3, 2, -7/6) © 
(3,1,-1/3) © (3,1,4/3) © (1,2,1/2) contains a weak 
triplet (A 3 ), doublet (A4), and singlet (A 5 ) leptoquarks 
we are interested in whereas the 5-dimensional Higgs rep- 
resentation, 5h, contains a singlet leptoquark only. So, 
these leptoquark states must be present in any renormal- 
izable theory based on SU(5). 

The most stringent constraints on leptoquark masses 
and their couplings to matter originate from limits on 
partial proton decay lifetimes. In that respect only A4 
is innocuous enough since it does not directly mediate 
proton decay. (It cannot couple to a quark-quark pair.) 
It is also practically impossible for it to be a part of 
the process that destabilizes proton through mixing with 
the Higgs doublet and some other state that couples to 
a quark-quark pair since (3,2,1/6) - the only suitable 
candidate - is not part of cither 5- or 45-dimcnsional 
Higgs representation (or 24-dimensional representation). 
It is thus phcnomcnologically possible for A4 to be light 
and have couplings to the matter fields of the form given 
in Eq. (5) in SU(5). In that case Y 2L = -2 1/2 [Y 4 \ - Y 4 %] 
and Y 2R = Y 4 \, . 

It is also possible to have A 3 that couples to the quark- 
lcpton pairs and no proton decay. In particular, the 10- 
5-45^ contraction yields a lepton-quark pair couplings 
with A3 of the form given in Eq. (2): Y 3 = Y 45 ,. On the 
other hand, the 10-10-45# contraction yields couplings 
of A 3 to a quark-quark pair only. Clearly, if only one 
of these two possible contractions is present there would 
not be a tree level proton decay due to A3. In the former 
case there would not be proton decay due to the mixing 
of A 3 with the Higgs doublet and some other states either 
since (3, 1, 2/3) - which would be a suitable candidate - 
is not part of cither 5- or 45-dimensional representation. 

Finally, A 5 could also be coupled to matter in a man- 
ner that renders proton stable contrary to the usual ex- 
pectation. Namely, the 10-10-45/f contraction yields 
couplings to a lepton-quark pair only. (This should be 
compared to the 10-10-5# contraction that generates 
both the lepton-quark and quark-quark type of couplings 
simultaneously for the singlet leptoquark in 5#.) This 
contraction yields the second term in Eq. (8): Ym — 
2 1 / 2 [Y 4 t 5 - Y 4 * 5 ]. The 10-5-45^ contraction, on the other 
hand, yields not only the second term in Eq. (5), i.e., 
Y\l = —2 1 / 2 Y£ 5 ,, but also the quark-quark couplings 
which would lead to proton instability. If only 10-10- 
45 h contraction is present proton could be stable and 
accordingly A5 could be light. Interestingly enough, it 
is possible to have a scenario in which there would not 
be any leptoquark induced proton decay. The necessary 
condition for this to happen would be the absence of the 
10-5-45tr contraction. 



